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Extragalactic background of high energy neutrinos arising from the interactions 
of cosmic ray protons with far-infrared extragalactic background radiation is calcu- 
lated. The main assumption is that the cosmic ray spectrum at energies higher than 
10 s GeV has extragalactic origin and, therefore, is proton dominated. All calcula- 
tions are performed with taking into account the possible cosmological evolution of 
,-£3 ■ extragalactic sources of cosmic ray protons as well as infrared-luminous galaxies. 
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1 Introduction 



The extragalactic backround of high energy neutrinos studied in this paper 
arises from collisions of high energy cosmic ray (CR) protons emitted by local 
sources (e.g., by active galactic nuclei) in extragalactic space, with extragalac- 
tic photons. We will not consider here the more exotic sources of extragalactic 
neutrino bachground (ENB) such as annihilations of topological defects (e.g., 
cosmic strings) or decays of hypothetical relic supermassive particles. So, ENB 
studied here began to be produced in relatively late epochs (from cosmolog- 
ical point of view) of expansion when the galaxies of different types and, in 
particular, the sources of CR's already existed. The first calculations of ENB 
were performed about twenty years ago [1,2,3]. Later, the detailed calcula- 
tions of differential spectra of ENB were done in works [4,5,6,7,8]. All these 
authors took into account the interactions of CR protons with relic photons 
(T S 2.7 K) only. 
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Evidently, other components of the extragalactic radiation background also 
must be taken into account in calculation of the extragalactic background of 
high energy neutrinos. The interval of photon wavelengths which is potentially 
most important from this point of view is (1 — 1000)/xm, i.e., the infrared 
region. At last fifteen years infrared astronomy developed very intensively 
(see, e.g., reviews of Franceschini [9] and Hauser and Dwek [10] ), and now 
the rather well-grounded calculation of ENB from interactions of cosmic rays 
with extragalactic infrared photons became possible [11,12]. 

Two main recent achievements of observational astronomy studying sources 
and fluxes of infrared photons in extragalactic space are sufficient for the 
present work: i) the extremely infrared-luminous galaxies have been discov- 
ered in the distant Universe (i.e., at substantial redshifts, z ~ 1) and their 
brightness distribution were measured and ii) the intense isotropic diffuse 
background radiation in the far-infrared, of extragalactic origin, have been 
revealed, and the spectral intensity of this background appeared to be of the 
same order as that of the optical extragalactic background. The conclusion 
from these observational facts is that some galaxies in the past should have 
been much more bright in the far-infrared than in the optical. It means that 
a large fraction of the energy emitted by distant galaxies should have been 
reprocessed by dust at long wavelengths, and that violent processes of star 
formation took place in the past in such massive systems. This evidence of the 
strong cosmological evolution of infrared extragalactic sources is very impor- 
tant for theoretical predictions of ENB as we show in this work. 

Technically, the problem of calculating the ENB from infrared photons is 
rather complicate because now the photon number density in extragalactic 
space, as a function of cosmological time, is not given (in the previous case of 
relic photons this number is definitely known for all moments of time) and, 
therefore, must be calculated separately. 

For the calculation of ENB one needs the following inputs. 

1. Photon energy spectra of infrared-luminous galaxies and the number density 
of these galaxies in space. Extensive information concerning these character- 
istics is available in literature [9,10]. For the concrete calculation of ENB we 
chose the model of Beichman and Helou [13] who derived the synthetic spectral 
energy distributions in the large diapason of values of the total luminosities, 
and the simple double power-law form of the local luminosity function (i.e., 
the present-day number density) of infrared galaxies suggested by Soifer et al 
[14] and based on the IRAS data [15]. 

2. The parameters characterizing the density and luminosity evolution of the 
infrared sources with look-back time. Although the fact of a strong time evo- 
lution of extragalactic infrared local sources is well grounded, it is reasonable, 
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however, at the present stage, to keep these parameters free. 

3. Extragalactic CR spectrum and intensity and the corresponding evolution 
parameters for the sources of CR 's. We do not know well enough the extra- 
galactic cosmic ray spectrum and, therefore, we are forced to use the crucial 
hypothesis about extragalactic origin of high energy CR's. Everywhere in our 
calculations we use the extragalactic model (the crossover energy ~ 10 17 eV) 
and normalize our theoretical CR spectrum on experimental CR data. 

Neutrino flux from interactions of high energy cosmic rays with infrared back- 
ground was recently calculated by Stanev [11]. Our work differs from that of 
[11] in several respects. The main difference is that in [11] there is no non- 
trivial cosmological evolution of infrared background: this evolution is assumed 
to be the same as of microwave background radiation (with z max = 8). Such an 
simplified approach requires minimum information from infrared astronomy: 
only the spectral intensity of the radiation at our epoch is needed, which can 
be, in principle, taken directly from experimental data. Besides, in [11] there 
is no separating of contributions to ENB from infrared and optical diapasons 
of the background while we calculate ENB from far-infrared part of the radi- 
ation background only (incidentally, it is well known that the time evolution 
of optical sources is much weaker in comparison with infrared ones). 

The paper is organized as follows. In Sec. 2 we derive approximate expression 
for the extragalactic CR spectrum for nonzero redshifts using the continuous 
energy loss approximation. Sec. 3 is dedicated to calculation of far infrared 
photon background and to obtaining the constraints on parameters of cosmo- 
logical evolution of this background. In Sec. 4 we calculate the extragalactic 
neutrino spectrum produced by interactions of cosmic ray protons with far 
infrared background, using the inputs from Sec. 2 and Sec. 3. The last section 
contains the result of the calculations. 



2 Extragalactic CR proton spectrum at different epochs 

To obtain approximate expressions for the spectra and intensities of CR pro- 
tons in extragalactic space we use the cosmological transport (kinetic) equa- 
tion without integral term, i.e., we work in the continuous energy loss approx- 
imation introduced, for these problems, by Berezinsky and Grigor'eva [16]. 
The validity of this approximation was studied in subsequent works, e.g., in 
the work of Yoshida and Teshima [7]. It was shown that this approximation 
works rather well and is not very good only for calculations of " local" cosmic 
ray spectra, i.e., cosmic ray spectra near the Earth from separate extragalactic 
sources (in particular, the spikes near the GZK cut off predicted in calculations 
using the continuous energy loss approximation are more sharply expressed 
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than in exact Monte-Carlo calculations). But in the case of diffuse cosmic 
ray spectra (i.e., the spectra integrated over redshift), the approximation of 
continuous losses is satisfactory: in particular, it predicts correctly the charac- 
teristic dip-bump structure of the cosmic ray spectrum between 10 19 and 10 20 
eV, at least in models with not very strong cosmological evolution of sources. 

The cosmological transport equation for extragalactic CR protons is in our 
approximation (see, e.g., the work of Blumental [17]) 

dn(E,z) d rn/rn . /rn ., 3n(E, z) n(E,z) , _. , /1X 

Here, n(E, z) is the differential number density of protons, having the di- 
mension (unit volume) -1 • (unit of energy) -1 , i.e., n(E,z)dE is the number of 
protons with a given redshift z in the energy interval dE near E. The function 
f3(E, z) is the change (loss) of proton energy in unit interval of z, 

„/ ,_, , dE E , . ,_, . dt 

The first term in r.h.s. of Eq.(2) takes into account adiabatic energy losses 
(those due to the cosmological expansion). The function b(E,z) describes the 
(non-adiabatic) proton energy losses per unit of time spent in extragalactic 
space filled by photons of the diffuse radiation background, and is defined by 
the expression 

t;\ lp ,z) = -b(E,z), (3) 

z 

where t~ 1 (7 p , z) is the cooling rate of protons with gamma-factor j p via jry — > 
nX and jry — > pe + e~ reactions at the cosmological epoch with redshift z. 
Evidently, the cooling rate, as a function of z, depends on the cosmologi- 
cal evolution of the radiation background. Throughout this work we use the 
approximation that the relic photon gas (i.e., only one component of the ra- 
diation background) is entirely responsible for the energy losses of cosmic ray 
protons in extragalactic space. In this approximation one has 

t~\ lpi z) = (l + zft-\ lp (l + z),0). (4) 

Using the text-book expression defining the energy loss coefficient b(E) , 

one obtains the connection of the function b(E,z), that enters Eq.(2), and 
b(E). 

b(E,z) = (l + z) 2 b(E(l + z)). (6) 

The function X(E,z) in the kinetic equation is connected with the mean free 
path X(E,z) of a cosmic ray particle for an interaction leading to its absorb- 
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tion. For simplicity, we neglect the corresponding term in the kinetic equation 
considering all proton interactions as non-catastrophic ones (in accordance 
with our approximation of continuous energy losses). Of course, some absorp- 
tion of cosmic ray protons in extragalactic space takes place, due to interac- 
tions with matter along the proton's path, and for taking it into account the 
absorbtion term must be kept. 

The function g(E, z) in r.h.s. of the kinetic equation describes the combined 
source of extragalactic cosmic rays. This source function can be written in the 
form 

g(E,z)=p(z) V (z)f(E)^. (7) 

Here, p(z) is the number density of local CR sources (e.g., AGNs) in the proper 
(physical) volume, 

p(z)=p (l + z) 3 , (8) 

r](z) is the activity of each local source (the integrated number of produced 
particles per second), 

V (z) = (l + z) m Vo 6(z max -z). (9) 

Writing Eq.(9) we assume that the cosmological evolution of cosmic ray sources 
can be parametrized by power law with the sharp cut-off at some epoch with 
redshift z max (m and z max are considered as parameters of a model of the 
combined source). At last, the function f(E) in Eq.(7) describes a form of the 
differential energy spectrum of the local source and has the simple normaliza- 
tion: 



/ f(E)dE 



1. (10) 

For the connection of differentials of cosmological time and redshift we use the 
simple formula 



dz H {1 + z) 



- 1/2 mi 



(assuming a flat Universe with matter and cosmological constant components). 
The matter density, in units of the critical density, Q m , is assumed to be equal 
0.3, and Q\ + Vt m = 1. 

The solution of the kinetic equation is given by the expression 

n(E,z)= J dxexp j / - ^^f ) dy\ ■ g(Z(x),x), (12) 
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and the function £ (y) is a solution of the equation 



d^y) _ ay) + (i + y)b(ay) ■ + y)) 



dy i + y H [n m (i + y y + n A } 1/2 

with the initial condition 

By differentiating the /3 function given by Eq.(2) one obtains 



dP(ay)-y) 



+ 



[i + yf 



day) i + y ' /J [fi m (i + y) 3 + fiA 

so, the final formula for the solution is 



(13) 



(14) 



(15) 



z 

n(E, z) = J dx ^ 



1-Kz 
1 + x 



X 



x cxp 



:i+y) s 



tfotO^l + ^ + Ov] 



1/2 



b\E) 



(16) 



^(£(x),x). 



Now, for the sake of the illustration, we consider the approximation which 
allows one to obtain the analytical solution of the kinetic equation expressed 
through the simple quadratures, which is rather good for approximate calcu- 
lations and estimates. 

Namely, we approximate the energy loss coefficient b(E) by the following sim- 
ple expressions: 



(17) 



b(E) = p E ; 3 • !0 9 GeV <E< !O u GeV, 

b(E) = (3' E ; E > 10 u GeV 

and b(E) = for E < 3 • 10 9 GeV. Constants (3 and (3' are: 

P' = 2 ■ IQ-Syear-- 1 ; (3 Q = 2.5 • lQ-^year' 1 . 



We assume, further, the power law for the injection spectrum from the extra- 
galactic local source, 



f(E) = J inj (E) 



^'^(GeV-i). 



E \E 



(18) 
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(E is the cut-off vafue for the extragalactic part of the spectrum of the local 
source; presumably, the particles with lower energies are somehow confined 
inside the source). 

At last, we put in Eq.(ll) Q m = 1, Q\ = 0. This approximation is not bad 
near z = 0; even at z — 3 it underestimates dt/dz on a factor which is smaller 
than 2. 

The resulting formulae for the cosmic ray proton spectrum appear to be dif- 
ferent in three energy regions (everywhere below proton energy is substituted 
in GeV): 



I. 
2. 
3. 



E < 3 • 10 9 /(1 +z) 
3 • 10 9 /(1 + z) < E < 10 n /(l + z) 
E > 10 u /(l + z). 



In all final expressions for the spectrum there is the following common factor: 
a(E, z , 1 , Eo) = ±.^.l^±(l + ,r« ■ (I.) " . (19) 



Ait Hq Eq 
In the region 3. one obtains (J(E, z) = -^n(E, z)) 



\EqJ 



J(E,z) = Q J dx(l +x) m -f- 7 exp{l78(l -7) [(1 + xf' 2 - (1 + zf' 2 ] } . 



(20) 



For the region 2. the formula for the spectrum is slightly more complicate: 

min(zmax,x b ) 



J(E,z) = n J dx(l + a;) m -i" 7 exp{2.2 



;i + xf 2 - (1 + zf 2 }} 

(21) 



and the value Xb is determined from the equation 

;i+a; fe ) 3 / 2 -(l + ^) 3 / 2 



ln^ 1 + fj=2.2 



E(l + x b ) 2 
In the region 1. there are two subregions: 

la. x bl = yj^-(l + z ) - 1 > z max ; E > E ; 

1 



(22) 



j(e,z) = n- 



m- \ -7 



i+z max r-^-(i+z) m ->- 



(23) 



lb. x bl < z„ 
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j(e,z) = n 



m — 5 — 7 



;i + x bl ) m -^ - (1 + z) m —^ 

(24) 



mi»(2mai,a:(, 2 ) 

+ f rfa;(l + a;) m -i- 7 exp{2.2(l-7) [{1 + x) 3/2 - (1 + x bl f 2 ]} 



and the value x b2 is determined from the equation similar to Eq.(22): 
10 n (l + z) 



In 



E{i + Xb2 y 



2.2 



(1 + x b2 f 2 - (1 + x bl f 2 



(25) 



Using expressions derived in this section one can calculate extragalactic cosmic 
ray proton spectra for different cosmological epochs. The parameters that 
are necessary for such a calculation are: m, z max , 7 and the product 
(PoVo)- Calculating J(E,0) = ■ n(E,0) with fixed values of m, z max , 7 
and normalizing it on the measured cosmic ray spectrum at energies 10 17 — 
10 19 GeV one can determine this product characterizing the combined cosmic 
ray source. As a result of this procedure we obtain the extragalactic cosmic 
ray proton spectrum for different epochs which can be used for the ENB 
calculation. 



3 Far infrared extragalactic background 



For a calculation of the extragalactic radiation background it is convenient to 
use the cosmological transport equation which is analogous to that used in the 
previous section. The function which must be found is the number density of 
infrared photons at different cosmological epochs, N IR (E 7 , z). The expression 
for the source function of the kinetic equation is 

9 IR {E^z) = J ^p(z,L)S IR (E„L)±- ■ f z (26) 

Here, the function S IR is the spectral luminosity ("spectral energy distribu- 
tion") of a luminous local source (which is in this case an infrared-luminous 
galaxy with a given total infrared luminosity L). The normalization of this 
function is 

J S IR (E J , L)dE 1 = L, (27) 

K 
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where, restricting oneself by the far infrared region, one takes 

E min = L 243-l(r 3 eV (A = 10 3 /im) ; E™ ax = 1, 243-KT 1 eV (A = 10pm). 

(28) 

The function p(z,L) in Eq.(26) is the number density of infrared-luminous 
galaxies with a given luminosity L in the proper (physical) volume, which is 
connected with the local luminosity function p(0, L) by the relation 

p(z,L)=p(0,^-^)(l + z)^. (29) 

where 7^ and 7/ are parameters determining the cosmological evolution of the 
luminosity of the local source, and the comoving density of these sources, 
respectively. 

For numerical calculation of the spectral energy distributions we use the model 
of Beichman and Helou [13]. These authors determined (theoretically) the syn- 
thetic spectral energy distributions F(E 1 , L) for infrared-bright galaxies with 
total far infrared luminosities in the interval (10 8 — 10 11 ) • L , taking into ac- 
count a number of components, the main of which are the "cirrus" component 
(infrared emission from a galactic diffuse dust heated by the interstellar radi- 
ation field) and a hot starburst component (originating from the violent star 
formation processes occuring inside dusty molecular clouds). The distributions 
F(E 1 , L) were calculated by Beichman and Helou in units of a spectral flux of 
energy (in Jy), for a reference distance 1 Mpc from the galaxy. The connection 
between F(E J ,L) and S IR (E 1 ,L) is 

S IR (E J , L) = 2 • 10 53 F(£ 7 , L), (30) 

and S IR (E^, L) has dimension s^ 1 (the luminosity per unit interval of energy). 

For the local luminosity function p(0, L) we use the analytic fit of Soifer et al 
[14] obtained from studies of IRAS infrared galaxies in the wavelength interval 
12 - 100/im, 

p(0,L) = 10- 2 - 07 -L- 2 - 31 , L>10°< 4 

(31) 

p(0,L) = 10- 2 - 73 -L-°- 65 , L<10 04 . 

In Eq.(31) the total luminosity L is measured in units of (10 10 L Q ) and p(0, L) 
is in (Mpc)~ 3 (In L)^ 1 . The luminosity function was modelled in the interval 

10~ L8 < L < 10 3 , 

and was set to zero outside this range. 
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For a derivation of the final expression for the number density of infrared 
photons in extragalactic space we use the formula (16) of sec. 2, putting there 
g = g IR and 6 = (i.e., we neglect energy losses of the infrared photons during 
their travelling in space). Substituting Eq.(26) into Eq.(16) one obtains 



Zmax _, , 2 



N IR (E„z)= J dx(j±^ J T p(x,L)S IR (^(x),L)x 

z 



(32) 



1 

x- 



The function £ 7 (a0 is a solution of the equation (compare with Eq.(13)) 

d^(x) _ £ 7 (x) 



dx 1 + x 



(33) 



with the initial condition 
and, therefore, 



£r(*) = E ^ (34) 

£,(*) = E,\±^. (35) 
1 + z 

The function N IR (E 1) z) is measured in cm~ 3 eV _1 if E 1 is substituted in (eV), 
p in cm _3 (lnL) _1 and S IR in s _1 . 

The expression (32) is the main analytic result of this section. It allows to 
calculate the extragalactic radiation background for any cosmological epoch 
with taking into account the time evolution of the luminosity as well as density 
of the local sources. In the particular case z = one has from Eq.(32) the 
number density of infrared photons for our epoch. Multiplying it on E 1 and 
on a factor ^ one obtains the well known formula [18] for the spectral flux of 
the extragalactic background radiation, 



I(E 1 ) = ±fL(R r (l + z),. 



dt 



dz 



dz. (36) 



Here, L(£ , 7 ,z) is the spectral luminosity density of all luminous sources in a 
comoving volume element at redshift z, 

L(£ 7 (l + *)) = J dz j d\nL Pc (z,L)S IR (E J (l + z)), (37) 

p c is the comoving luminosity function 

^• i) = frrf' (38) 
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The function E 1 I(E 1 ) is the "spectral intensity" of the radiation background 
(having dimension (energy-cm _2 -sec _1 -ster -1 )), which is commonly used for a 
presentation of data of background measurements. 

As is mentioned in the Introduction, a distinctive property of infrared-luminous 
galaxies is their very high rates of cosmological evolution, i.e., evolution with 
redshift z. These rates sizeably exceed those measured for galaxies at other 
wavelengths. In particular, the comoving luminosity density from luminous 
infrared galaxies was at z ~ 1 on a factor about 10 (or even more) larger than 
in the local Universe. Correspondingly, the theoretical results which can be 
obtained for the extragalactic infrared radiation background are very sensi- 
tive to chosen values of the evolution parameters 7/, 7^. From data on number 
counts of infrared galaxies at 15 and 24 \im one can conclude that there is 
some evidence of a strong density evolution as well as strong luminosity evo- 
lution (7; ~ 7^ ~ 3), although one cannot still exclude the possibility of an 
almost pure density evolution (7^ 3> 7;) 

The detailed analysis of constraints on the evolution parameters following 
from data is given in the review of Franceschini [9]. It is argued in [9] that 
the assumption of the very high evolution rates for the population of infrared 
sources (up to 7^ ~ 8 or more) is not in contradiction with data. The strong 
density evolution seems to be very natural, at least for the population of star- 
burst galaxies (those dominating at the high luminosity end of the luminosity 
function) because the starburst phenomenon is caused, most likely, by interac- 
tions and merging of galaxies (and the probability of such interactions scales 
as square of the number density of galaxies in the physical volume). 

In the concrete numerical calculation we used the following walues of the 
evolution parameters 

7 d = 4,6,8 ; 7; = 0, ; z max = 1, (39) 

i.e., we assume that there is no luminosity evolution and that there is the sharp 
cut-off after z max . On Fig.l we show the results of our far infrared background 
calculation for our epoch, together with experimental data, taken from [9]. We 
remind the reader, that we restrict ourselves in this work by considering only 
the far-infrared part of the background (10 — 1000/xm). It is seen from the 
figure that the choice 7^ = 6 looks as the most appropriate. The slight shift of 
the maximum of the curves aside from the data is entirely due to the features 
of theoretical spectral energy distributions used in the calculation. 

One should note that our choice of the parameters 7^, z max is the simplest 
one. In the more realistic situation there is no sharp cut-off of the evolu- 
tion: it is most probable, from physical point of view, that the evolution rate 
changes gradually. Any possible variants can be considered modifying, triv- 
ially, Eq.(32). 
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Fig. 1. Extragalactic muon neutrino spectrum from interactions of cosmic rays with 
far infrared background radiation, for two values of density evolution parameter 7^. 



4 Extragalactic spectra of high energy neutrinos 



For simplicity we use in this work one-pion approximation, i.e., we assume 
that in jry-reaction only two particles (neutron and pion) are produced. In 
this case, as is well known (see, e.g., Hill and Schramm (1985)) the pions have 
approximately isotropic distribution in the center of mass system and, as a 
consequence, the step-like energy spectra in the observer system. 

The main photoproduction reaction, 

PCR + linf raved -> K + + TO, (40) 

and subsequent decays 

7r + — > fi + + Vp, fi + -> e + + v e + i> M (41) 

lead to production of v^v^ u e . In the present work we will be interested only 
in (z/ M + z> M )-flux. Therefore, we add together the neutrino from pion and muon 
decays. Neutron produced in the photoproduction process will decay giving 
proton, long before the next interaction with photons of the background. 

The number of proton-photon collisions per second and per units of photon 
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energy and cos 6* 7 is equal to (for a given z) 

dn coll c lti 



= -a(s)(l + cos O^N^iE^z). (42) 



dtdE^d cos # 7 2 

Here, er(s) is the total cross section of the reaction (40). Using the connection 
between £ 7 and E l ® b , 

E\ ab = £ 7 (l + cos0 7 ) - 7p , (43) 

one obtains 



f/ "'"'' - ^-a(Elf)N IR (—J%lL-—, z ) . (44) 



dtdE^d cos 9^ 2^ v K 7 7 \7 P (1 + cos 7 ) ' ' 

Evidently, n C0 22 is a function of -E 7 afc , cos# 7 and z. Differential neutrino 
production spectrum from the pion and muon decays is given by the following 
integral: 



= iJdE p f iE? J i cos ft, ME £' m - M E r , E?, E l ,)J(E pl z) l 

(45) 

where J(E p , z) is the extragalactic CR proton spectrum calculated in sec. 2 and 
ju{E p , E 1 ^, E v ) is the neutrino spectrum per one p 7 collision. The production 
spectrum is measured in (cm^sec^GeV^ 1 ), and, as it should be, it is the 
number of neutrinos produced in unit volume per second and per unit interval 
of energy. So, this value has the same sense as the product p(z)r)(z)f(E) in 
the case of CR protons (see Eq.(7)). 

The last step is the calculation of the neutrino extragalactic spectrum (or 
number density), i.e., the ENB, by integration the NP rod over all redshifts. 
The source function of the neutrino transport equation is 

g„(E„,z) = Nr d (E„,z)^ z (46) 
and the final result for the neutrino spectrum is 



UEu) = ^ J (1 + z)2 g„(E„(l + z),z)dz = 

(47) 

_ l/ -Nr d (E u (l + z),z)dz. 



-max 



c 
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5 Results and discussions 



As is noted in the Introduction, we used in this work the basic assumption 
that ultra high energy cosmic rays have extragalactic origin. It means that 
at E > E , where E is some crossover energy, all cosmic rays detected near 
the Earth came from extragalactic space where they had interactions with 
extragalactic radiation background and, moreover, they are mostly protons. 
This hypothesis allows us to normalize our theoretical cosmic ray spectrum 
obtained in sec2 on data of cosmic ray measurements at energies > E Q . 

We did not attempt in this work to sew accurately the galactic and extra- 
galactic spectra near crossover energy, it is the separate and delicate problem. 
It is enough for us, at this stage, that the characteristic features of the exper- 
imental cosmic ray spectrum at 10 18 — 10 19 eV, i.e., the dip-bump structure 
and the beginning of the GZK cut-off, can be well described in our model, 
by the proper choice of the parameters 7,m, z max . In general, the possibility 
of a describing of the cosmic ray spectrum at E ~ 10 18 — 10 19 eV by models 
with extragalactic origin of high energy cosmic rays was proved in a series of 
papers by Berezinsky with coauthors [16,19]. 

We chose for the crossover energy the value E = 10 17 eV. The final results do 
not depend very much on the concrete choice of E if E Q is much lower than 
the threshold of the photoproduction on infrared photons, 

~ mpmn ~ l0 19 eV. (48) 

1 far—infr 



The experimental form of the cosmic ray spectrum at 10 18 — 10 19 eV is well 
described with the following set of parameters: 

7 = 2.5 , m = 3.5 , z max = 5. (49) 

We see that the value of z max characterizing the evolution of cosmic ray 
sources is much higher than the corresponding value chosen in sec. 3 for infrared 
sources. It is not the contradiction because the physical nature of sources is 
clearly different in these two cases: the contributions of AGNs to the infrared 
background is rather small (see, e.g., [20]). 

Normalizing the theoretical cosmic ray spectrum at z = with the chosen set 
of parameters on cosmic ray data one obtains 

Po r] 10- 41 (cm _3 sec _1 ). (50) 

This value corresponds to the total luminosity required for the production of 
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Fig. 2. Far infrared radiation background calculated using the model of Beichman 
and Helou, for three values of the density evolution parameter 7^. 



cosmic rays with E > E Q 



4.5 • 10 45 (ergMpc^Yr- 1 ). 



(51) 



If, e.g., there is one cosmic ray source per volume 10 6 (Mpc) 3 , the activity of 
each such source is about 10 44 erg/s (for E > 10 17 eV). 

The resulting curves of the muon neutrino extragalactic background from in- 
teractions of cosmic rays with far- infrared photons of the radiation background 
are shown, for two chosen values of the parameter 7^, on Fig. 2. One can see a 
strong dependence of the result on the assumed evolution law of the infrared 
background. The order of magnitude of ENB is qualitatively the same as in 
the Stanev's work [11]. In spite of the fact that the density of infrared pho- 
tons in extragalactic space is much smaller (~ 1.5 photons/cm 3 ) than that of 
relic microwave photons, the neutrino background appears to be not so small, 
due to the lower threshold for photoproduction and, last but not least, due to 
much stronger time evolution of infrared background in comparison with that 
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of relic photons. The predicted neutrino fluxes are comparable, more or less, 
with the neutrino fluxes from other extragalactic sources at energy region near 
10 17 eV (7-ray bursts, topological defects, etc) and deserve further theoretical 
and, in future, experimental studies. 
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